Spin-orbit coupled particle in a spin bath 
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We consider a spin-orbit coupled particle confined in a quantum dot in a bath of impurity spins. 
We investigate the consequences of spin-orbit coupling on the interactions that the particle mediates 
in the spin bath. We show that in the presence of spin-orbit coupling, the impurity-impurity 
interactions are no longer spin-conserving. We quantify the degree of this symmetry breaking and 
show how it relates to the spin-orbit coupling strength. We identify several ways how the impurity 
ensemble can in this way relax its spin by coupling to phonons. A typical resulting relaxation rate 
for a self- assembled Mn-doped ZnTe quantum dot populated by a hole is 1 /is. We also show that 
decoherence arising from nuclear spins in lateral quantum dots is still removable by a spin echo 
protocol, even if the confined electron is spin-orbit coupled. 

PACS numbers: 75.30.Et, 76.60.-k, 71.38.-k, 73.21. La, 33.35.+r 
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I. INTRODUCTION 

A singly occupied quantum dot is the state of the art of 
a controllable quantum system in a semiconductor Co- 
herent manipulation of the particle spin has been demon- 
strated in lateral dots, where top gates allow for as- 
tonishing degree of control by electric fields^— and in 
self- assembled dots, where a weaker control over the dot 
shape and position is compensated by the speed of the 
optical manipulationi 7 - In both of these major groups, 
there are two main spin-dependent interactions of the 
confined particle and the semiconductor environment: 
spin-orbit coupling embedded in the band structure, and 
spin impurities, which are either nuclear spins, or mag- 
netic atoms 

A particle couples to an impurity spin dominantly 
through an exchange interaction, which conserves the 
total spin of the pairJ^ This way the electron spin 
in a lateral GaAs dot will decohere within 10 ns due 
to the presence of nuclei*^— Typically, such decoher- 
ence is considered a nuisance that can be partially re- 
moved by spin echo techniques prolonging the coherence 
to hundreds microseconds J£r— Whether that decoher- 
ence time can be extended further, e.g. by polarizing 
the impurities ) 21 ' 22 is not clear, as the experimentally 
achieved degree of polarization has been so far insuffi- 
cient, despite great effortsi 2 ^— On the other hand, a 
strong particle-impurity interaction is desired in magnet- 
ically doped quantum dots.— ~— Here the confined par- 
ticle is central for both supporting energetically, and as- 
sisting in creation, the desired magnetic order of the im- 
purity ensemble. In fact, similar magnetic ordering can 
be traced to the studies of magnetic polarons in bulk 
semiconductors, for over fifty years 4^ The formation of 
a magnetic polaron can be viewed as a cloud of local- 
ized impurity spins, aligned through exchange interaction 
with a confined carrier spin42~— 

The conservation of the spin by the impurity-particle 
interaction is a crucial property. For example, the spin 



relaxation of the impurity ensemble is impossible with 
only spin-conserving interactions at hand. This moti- 
vates us to consider possibilities to break this symmetry. 
The first and obvious candidate is the spin-orbit cou- 
pling (SOC)i&ii Despite being weak on the scale of the 
particle orbital energies, it dominates the relaxation of 
the particle spin in electronic dots, as is well known^S 
because it is the dominant spin-non-conserving interac- 
tion. The questions we pose and answer in this work are: 
assuming the particle is weakly spin-orbit coupled, how 
strong are the effective spin-non-conserving interactions 
which appear in the impurity ensemble and what is their 
form? Is the induced particle decoherence still removable 
by spin echo? Is the particle efficient in inducing impu- 
rity ensemble spin relaxation, thereby limiting the achiev- 
able degree of the dynamical nuclear spin polarization^ 7 -? 
Can the magnetic order be created through the spin-non- 
conserving particle mediated interactions — that is, is 
this a relevant magnetic polaron creation channel? 

To address these questions we develop here a frame- 
work allowing us to treat different particles and impu- 
rity spins in a unified manner. We apply our method 
to two specific systems: a lateral quantum dot in GaAs 
occupied by a conduction electron with nuclear spins 
of constituent atoms as the spin impurities, and a self- 
assembled ZnTe quantum dot occupied by a heavy hole 
doped with Mn atoms as the spin impurities (readily in- 
corporated as Mn is isovalent with group-II atom Zn). 
Both of these systems are quasi-two dimensional, the par- 
ticle spin-orbit coupling is weak compared to the particle 
orbital level spacing and the particle-impurity interac- 
tion is weak compared to the particle orbital and spin 
level spacingS i 2 ' 32 ' 48 ' 49 As it is known, 50 in this regime 
one can derive an effective Hamiltonian for the impu- 
rity ensemble only, in which the particle does not ap- 
pear explicitly. This can be done including the particle- 
bath interaction perturbatively in the lowest order, see 
the scheme in Fig. 1. Our contribution is in show- 
ing how the procedure generalizes to a spin-orbit cou- 
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pled particle. In addition, we use the resulting Hamil- 
tonian for the calculation of the spin relaxation of the 
impurities which is phonon-assisted (required to dissi- 
pate energy) and particle-mediated (required to dissipate 
spin). We come up with (and evaluate the correspond- 
ing rates for) five possible mechanisms how the spin flips 
can proceed: shifts of the particle by the phonon electric 
field (Sec. IV. B), position shifts of the impurity atoms 
(Sec. IV.B), relative shifts of bulk bands (Sec' IV.C), 
rcnormalization of the spin-orbit interactions due to band 
shifts (Sec. IV. D) and spin-orbit interactions arising from 
the phonon electric field (Sec. IV. E). 



Our main findings are the following: 1. The spin- non- 
conserving interactions couplings are proportional to the 
spin-conserving ones multiplied by some power of small 
parameter(s) which quantify the spin-orbit interaction. 
For the electron, the small parameter is the dot dimen- 
sion divided by the spin orbit length and the proportion- 
ality is linear. For the hole, the small parameters are the 
amplitudes of the light hole admixtures into the heavy 
hole states. The proportionality differs (from linear to 
quadratic) depending on which hole excited state medi- 
ates the interaction. The interaction form is given in 
Eq. (j43|) . our main result. 2. For the electron, the addi- 
tional decoherence is removed by the spin echo, while for 
holes only a partial removal is possible. The latter is be- 
cause, unlike for electrons, the spin-non-conserving cou- 
pling is mediated rather efficiently through higher excited 
states. 3. The piezoelectric acoustic phonons are most ef- 
ficient in relaxing the impurity spin. The resulting relax- 
ation time is unobservably long for nuclear spins, while 
the hole- induced Mn spin relaxation time of 1 fj,s is typical 
for a 10 nm self- assembled quantum dot, where experi- 
mentally measured times for the polaron creation range 
from nano to picoseconds! 26 ' 31 ' 32 From this we conclude 
that the interplay of spin-orbit coupling and phonons 
does not govern the dynamics of magnetic polaron forma- 
tion at moderate Mn densities (few percents) , but rather 
represents the spin-lattice relaxation timescale, similarly 
as is the case in quantum wellsi 41 ' 51 Apart from a very 
low Mn doping, the analytical formulas presented in this 
work allow us to identify additional regimes where the 
particle mediated spin relaxation could be relevant for 
the polaron creation: an example is a hole located at a 
charged impurity. 



The article is organized as follows: In Sec. II we in- 
troduce the description of the particle focusing on the 
spin-orbit coupling. In Sec. Ill we specify the particle- 
impurity interaction, define its important characteristics 
and derive the effective Hamiltonian for the impurity 
ensemble. In Sec. IV we calculate the spin relaxation 
rates for the impurity ensemble, after which we conclude. 
We put numerous technical details into Appendices, with 
which the text is self-contained. 
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FIG. 1: (Color online) Effective interaction between impuri- 
ties (encircled in red/gray) mediated by a confined particle 
(red/gray lines), a) Electron is excited from the ground state 
of spin +1/2 into the closest spin -1/2 state (up by the Zeeman 
energy A z ) upon flipping one of the impurities and deexcited 
back upon flipping another one. b) Hole spectrum is more 
complicated, comprising heavy hole (hh) and light hole (lh)- 
like states, the latter displaced by light-heavy hole splitting 
A th . 



II. QUANTUM DOT STATES AND SPIN-ORBIT 
INTERACTION 



Electron states 



In the single band effective mass approximation, that 
we adopt, the Hamiltonian of a quantum dot electron is 
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ffdot = ^-+V(r) + ^+V z (z)+g^ B B-3 + H so . (1) 
2m 2m 



The underlying bandstructure is taken into account as a 
renormalization of the mass m and the g factor in the 
electron kinetic and Zeeman energies, respectively. The 
latter couples the external magnetic field B to the elec- 
tron spin through the vector of Pauli matrices <x = 2 J. 
We will below assume a sizeable (above 100 mT) ex- 
ternal magnetic field in the electronic case, which is 
typical in experiments to allow for the electron spin 
measuremen t 46 ' 52 and to slow down the impurity dynam- 
ics and the resulting decoherence i^ 6 - On the other hand, 
we neglect the orbital effects of this field which is justi- 
fied if the confinement length is much smaller than the 
magnetic length lg = ^2h/eB, where — e is the electron 
charge. If the field is strong (above 1-2 Tesla), the orbital 
effects important here arc fully incorporated by a rcnor 
malization of the confinement length l~ 4 —> l~ A + I 
The field is applied along So, a unit vector. 

The quantum dot is defined by the confinement poten- 
tial V + V z , which we separated into the in- plane and 
perpendicular components. The corresponding in-plane 
and perpendicular position and momentum components 
read R = {r, z} and P = {p,p z }, respectively. When- 
ever we below need an explicit form of the wavefunction, 
we assume, for convenience, a parabolic in-plane and a 
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hard-wall perpendicular confinements 

V(r) = ^ . \ m .\\ (2a) 

V z {z) = (°' °< z<w ' ( 2b ) 
' 1 oo, otherwise. v y 

The confinement lengths / and w characterize a typical 
extent of the wavefunction in the lateral, and perpendic- 
ular directions, respectively. Alternatively to the length, 
the confinement energy hw can be used. However, we 
stress that our results below do not rely on the specific 
confinement form in any way, as long as the dot is quasi 
two dimensional, a condition which for the adopted ex- 
ample reads I 3> w. Typical values for lateral quantum 
dots in GaAs are /=30 nm, w=8 nm. 

The last term in Eq. (JlJ is the spin-orbit interaction^ 

H so = -^-{-VxPx + VyPy) + -j {&xPy - CJ y p x ), (3) 
ZTOij ZTOibr 

comprising the Dresselhaus term, which arises in zinc- 
blende structures grown along [001] axis and the 
Bychkov-Rashba term, which is a consequence of the 
strong perpendicular confinement. The interactions are 
parameterized by the spin-orbit lengths Z so <E {ld,lbv}, 
typically a few microns in GaAs heterostructures. 

Assume first that the spin-orbit coupling is absent. To 
be able to treat the electron and the hole (each referred 
to as the particle) on the same footing, we introduce the 
following notation 

\* P ) = \J)®\$ J a ), (zeroSOC). (4) 

The complete particle wavefunction, for which we use 
the Greek letter "J, is a two (electron) or four component 
(hole) spatially dependent spinor. Its label p — {J, a} 
indicates that the wavefunction is separable into a (posi- 
tion independent) spinor | J) and a scalar position depen- 
dent complex amplitude |<i>). The former is labelled by 
the particle angular momentum in units of h, J — ±1/2 
(electron), J = ±3/2, ±1/2 (hole; alternatively, we use hh 
for 3/2 and Ih for 1/2 labels). The set of orbital quan- 
tum numbers a depends on the confinement potential. 
For the choice in Eqs. ([2]), it is a set of three numbers 
a = {nm, fc}, with n the main and m the orbital quan- 
tum number (to = —to) of a Fock-Darwin state, and k 
the label of the subband in the perpendicular hard-wall 
confinement. Finally, for the particle ground state we 
omit the index a, or use G = {J, 00, 0} in place of p. The 
electron ground state is thus denoted by 

|*i/ 2 > = |1/2)®|#g), (5) 

where the direction of the spin up spinor |l/2) is set by 
the external field along So- 

Let us now consider the spin-orbit coupling. It turns 
out that for electrons the spin-orbit effects on the wave- 
function can be in the leading order written a o 53 i 54 

\9 p ) = U\J)<8>\$t). (6) 



Here U is a unitary 2 by 2 matrix of a spinor rotation, 

[/(r)=exp[-in so (r)-J], (7) 
parameterized by an in-plane position dependent vector 

Mr) = -(f-f,f-f,0). (8) 

A weak spin-orbit coupling allows us to label states with 
the same quantum numbers as for no spin-orbit case, as 
there is a clear one to one correspondence. The enormous 
simplification, that the unitary matrix in Eq. ([5]) does not 
depend on the quantum numbers p, is due to the special 
form of the spin-orbit coupling in Eq. @ , what has sev- 
eral interesting consequences To calculate the spin 
relaxation or spin-orbit induced energy shifts, on which 
U has no effects, one has to go beyond the leading order 
given in Eq. ([6]). 59 However, we will see that here U will 
result in spin-non-conserving interactions, and it is thus 
enough to consider the leading order. For the same rea- 
son we neglected the cubic Dresselhaus term in Eq. 
what is here, unlike usually^ an excellent approxima- 
tion. 



B. Hole states 

For holes we restrict to the four dimensional subspace 
of the light and heavy hole subbands. Neglecting the 
spin-orbit coupling, they correspond to the angular mo- 
mentum states J = ±3/2, and J = ±1/2, respectively. 
We use the confinement potential in Eqs. ([2]), setting the 
confinement energy in the heavy hole subband to 20 meV, 
which gives I ~ 4 nm and setting the light-heavy hole 
splitting to Aih = 100 meV, which gives w ~ 2 nm. The 
atomic spin-orbit coupling manifests itself as the orbital 
splitting of the light and heavy holes from the spin-orbit 
split-off band (which is energetically far from the states 
considered in this paper), and as a coupling of the light 
and heavy holes at finite momenta. The latter effect we 
refer to as the hole spin-orbit coupling — we do not con- 
sider higher order effects, which give rise to spin-orbit 
interactions similar in form to the electronic Dresselhaus 
and Rashba terms. Within this model, we derive the 
spin-orbit coupled wavefunctions from the corresponding 
4 by 4 sector of the Kohn-Luttinger Hamiltonian in Ap- 
pendix [A] and get the hole ground state as 

l*3/ 2 > =|3/2) ® |$ t ) + AJl/2) ® 
+ A |-1/2)®K )0 ). 

We use the notation explained below Eq. The Fock- 
Darwin states in the heavy hole (hh) and light hole (Ih) 
subband differ, due to different effective masses. The 
key quantities are the scalars A, which quantify the light- 
heavy hole mixing. The spin-non-conserving interactions, 
as well as the resulting spin relaxation rates, will scale 
with these scalars. For our parameters, which we list for 
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convenience in Appendix [El we get Ao ~ Ai 0.05. In 
addition to the ground state, we will need also the lowest 
excited state in the heavy hole subband, which is the 
time reversed copy of Eq. 

|*_ 3 /2> -I - 3/2) ® |$g* ) - A x | - 1/2) ® 1$"^) 
+ Ao|l/2)<8|${* )0 ), 

and also the lowest state in the light hole subband 
|* 1/3 ) =11/2) ® |< . ) + A;|3/2) ® |^ T \) 

11 

-A&|-3/2)®|$& )0 ), 

which will be surprisingly effective in inducing the spin- 
non-conserving coupling among impurities, as we will see. 
We make a few notes here: First, in the spherical approx- 
imation that we adopt, the Kohn-Luttingcr Hamiltonian 
conserves the angular momentum, so that all components 
in each of Eqs. ©-(HU have the same value of J + m. 
Second, the mixing is stronger in the light hole subband, 
X[ sa 0.15 and Aq 0.11. This is because the admixing 
states are closer in energy — the Fock-Darwin excitation 
energies add to and subtract from the light-heavy hole 
splitting in the heavy and light hole case, respectively, 
as evidenced by Eqs. (jA9l) and (|A11|) . Third, we will be 
interested in the case of zero external magnetic field for 
holes (we give the Zeeman interaction in Appendix [A] for 
completeness.) Unlike for electrons, such a field is here 
not required to split the two states in Eqs. © and (ITU1) . 
as the splitting arises due to the spin impurities. As we 
will see, this splitting will be of the order of few meV. 
Compared to this, the hole Zeeman energy is negligible 
up to fields of several Tesla. In addition, the external field 
suppresses an interesting feedback between the particle 
and impurities Finally, we note that one could relate 
the first-order and unperturbed hole states analogously 
to the electron case introducing a unitary transformation 
U, whose matrix elements are the coefficients appearing 
in Eqs. (l9|- (fTT]) . However, since here the transformation 
does not have any appealing form similar to the one in 
Eq. ([7]), we do not explicitly construct the matrix U for 
holes. 



III. EFFECTIVE HAMILTONIAN 

In this section we introduce the particle-impurity ex- 
change interaction, while providing a unified description 
for both type of particles and impurities. This interac- 
tion manifests itself as the Knight field acting on the 
impurities and the Overhauser field acting on the parti- 
cle. (The fields are defined as the expectation value of 
the exchange interaction in the state of the correspond- 
ing subsystem). Historically, the terminology was ini- 
tially applied to nuclei and here we also use it for Mn 
spins. With the help of these fields, we define the unper- 
turbed basis of the particle-impurity system, for which 
we derive the effective interaction Hamiltonian treating 



the non-diagonal exchange terms perturbatively. Finally, 
we define the spin-conserving versus spin-non-conserving 
interaction terms and analyze their relative strength for 
the electron and hole case. 

Our strategy can be viewed also in the following al- 
ternative way. To derive the spin-orbit coupling effects 
on the effective impurity interactions, we proceed in two 
steps: First we unitarily transform the particle basis 
to remove the spin-orbit coupling in the lowest order. 
The spin-orbit coupled basis transformation rcnormalizes 
the particle-impurity exchange interaction and breaks its 
spin-rotational symmetry. Then we integrate out the par- 
ticle degrees of freedom by a second unitary transforma- 
tion, using the Lowdin (equivalently here, the Schrieffer- 
Wolff) transformation, which leaves us with effective in- 
teractions concerning impurities only. 

A. Particle-impurity interaction 

The particle interacts with impurities by the Fermi 
contact interactio n 47 ^ 2 

H F =J2 H F =-^/3^(R-Rn)J-I". (12) 

n n 

Here n = 1,2... labels the impurities located at posi- 
tions R„ = (r n , z n ) with corresponding spin operators 
I™ in units of h. The impurities have spin I and den- 
sity 1/vq. For the electron- nuclear spins case the impu- 
rity density is in GaAs the atomic density. For holes-Mn 
spins we assume XMn is the fraction of cations replaced 
by Mn atoms, typically iMn = 1%- For impurities with 
different magnetic moments (as for nuclei of different el- 
ements) the coupling /3 should have the index n, but we 
will not consider this minor complication. Even though 
the wavefunction extends to infinity, it makes sense to 
consider the number of impurities with which the par- 
ticle interacts appreciably as N — V/vq, where the dot 
volume V is defined byi^ 

1/V= J d 3 R|$ G (R)| 4 . (13) 

The maximal value of the Hamiltonian Hp, if all impu- 
rities are aligned with the particle, is 

E = -J2P\®G(Rn)\ 2 JI = -l3JI/v . (14) 

n 

The impurity Zeeman energy is 

H nZ = ^2 3nMimpI" ' B, (15) 

n 

with g n the impurity g-factor. For GaAs lateral dots, 
we consider nuclear spins as impurities, with / = 3/2 
and /Ltimp the nuclear magneton //jy. For ZnTe dots, the 
impurities are intentionally doped Mn atoms, with I = 
5/2 and /Xi mp the Bohr magneton hb- 
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B. Knight field 

Assume the particle sits in the ground state G. In the 
lowest order of the the particle-impurity interaction, a 
particular impurity spin couples to a local field, called the 
Knight field. We define it in units of energy by writing 

K n -I n = {^ G \H^ G ), (16) 

from which, using Eq. (fT2]l , we get 

K" = - J 8<* G |«(R-R n )J|* G ). (17) 

Using Eq. (|6]) the Knight field of an electron is 

K" = -/3|$ G (R n )| 2 (l/2|C/t(R„)J[/(R„)|l/2) 
= -/3(l/2)|$ G (R„)| 2 (s(R„)). 

It points along the direction of the electron spin at the 
position of the n-th impurity, introduced as a unit vector 
(s(R n )) = i?[/( Rn )[so]. The operator Rjj is defined such 
that it performs the same rotation on vectors, as U does 
on spinors. The explicit form of R is the one in Eq. ([?]). 
if generators of rotations in three dimensions are used, 
(Jk)im — — iefeim- As is apparent from Eq. (fTS|) . evalu- 
ating the Knight field with perturbed electron wavefunc- 
tions is equivalent to evaluating a unitarily transformed 
interaction Hp — > WHpU with unperturbed electron 
wavefunctions. 

We get the Knight field of the hole as (see Appendix 
0, 

[K2,Kf,K2] = -flAxRe/n, AiIm/„, 3|*$ (R„)| 2 /2], 

(19) 

where we abbreviated /„ = \/3$o^,i(R«)*oo,ci(^n) an d 
neglected contributions of higher orders in A. Rather 
than the exact form, we note that without the spin- 
orbit coupling, the direction of the Knight field is fixed 
globally (along the external magnetic field for electrons, 
So = B/£? and along the z axis — the spin direction of 
heavy holes — for holes, So = z). The spin-orbit inter- 
action deflects the Knight field in a position dependent 
way, the deflection being in the leading order linear in the 
small parameter characterizing the spin-orbit interaction. 
In this respect, Eq. (IT5|) and (fTU|) are the same. 

C. Basis 

The total field aligning the impurity spin is the sum of 
the Knight field and the external field 

B" =K n + g nfMmp B. (20) 

The typical energy scale of the Knight field of an electron 
in a lateral dot is tens of peV, which corresponds to the 
impurity in an external field of 1 mT. For a hole in a self- 
assembled dot the Knight field is of the order of 100 /J,eV, 
corresponding to the external field of 0.3 T. Based on 



this, in the following analysis we mostly consider typical 
situations, in which the total field is dominated by the 
external field for nuclear spins (electronic case), and the 
Knight field for Mn spins (hole case). 

We now introduce for each impurity a rotated (primed) 
coordinate system, in which the unit vector z' is along 
the total field B n . Formally, the rotation is performed by 
operator i?B>» defined by the relation between the unit 
vectors, 

z' = i?B-N- (21) 

The orientation of the in-plane axes x' ', y' in the plane 
perpendicular to z' is arbitrary, and we denote r'± = 
x' ± iy'. We define the impurity ensemble basis states as 
tensor products of states with a definite spin projection 
along the locally rotated axis z', 

\l) = \ll,)®\ll)®...®\I%). (22) 

The spin projections take discrete values, I™, € {1,1 — 
1, . . . , — /}. We use 1 as the collective index of the im- 
purities. With this, we are now ready to introduce the 
complete system basis, as spanning the states 

\%) ® \T) = \%) ® |ji ) ® |J|,) ® ... ® (23) 

with the corresponding energies 

E PtI = E P + Y / % = E p + B ni z' . ( 24 ) 

n n 

comprising the particle energy and the Zeeman energies 
of impurities in the corresponding total fields. 

D. The substantial gap assumption — Overhauser 
field 

In addition to the Knight field, another consequence 
of the particle-impurity interaction from Eq. (|12p is the 
effective field experienced by the particle spin, known as 
the Overhauser fielcU O. To express this field again in 
the units of energy, it is helpful to consider the matrix 
elements of the particle-impurity interaction within the 
subspace of the lowest two electron states J, J' € S = 
{1/2,-1/2}, 

(*j|-/35>( R - R «) J - in l*->'} 

™ (25) 
= -/3^|$ G (R n )| 2 i? t/ t[F]-(J|J|J'). 

n 

We introduce the field O as 

Hf\s = O • J, (26) 

where the subscript S refers to the subspace comprising a 
pair of time-reversed particle states and the Overhauser 
field depends on the choice of S. To quantify the Over- 
hauser field, we give up on trying to track the microscopic 
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state of the impurities and instead introduce the averag- 
ing (denoted by an overline) over impurity ensembles 

7" = 0, W = S nm Sa b I(I + l)/3, (27) 

which characterizes unpolarized and isotropic ensembles. 
Nuclear spins, unless intentionally polarized in dynami- 
cal nuclear polarization schemes , 22 i 63 i 64 usually well fulfill 
this. If not polarized by an external field, Mn spins can 
be considered random initially, before the particle enters 
the dot and the polarization starts to build up. 

Equation (|27l) gives a zero Overhauser field on average, 
but with a finite dispersion, quantifying a typical value. 
For electrons, we get the well known result ) 13 ' 65 

o^ = (3 2 Yl l<MR„)l 2 l<MR m )l : %t [I"] • R uL [i m ] 

nm 

= p 2 J2 l<MR„)l 4 /(/ + 1) = J{i + l)(P/v ) 2 /N, 

n 

(28) 

stating that the typical value of the Overhauser field is 
inversely proportional to the square root of the number of 
impurity spins within the dot. The spin-orbit coupling, 
equivalent to a position dependent spin coordinate frame 
rotation, does not influence the result at all, as Eq. ([27]) 
assumes isotropic non-interacting impurities. For our pa- 
rameters, the typical Overhauser field value is 0.3 ^eV, 
which corresponds to external field of 20 mT. The energy 
splitting of the electron spin opposite states is therefore 
for our case dominated by the Zeeman, rather than the 
Overhauser, field. 

For holes, we will not write the Overhauser field ex- 
plicitly as a vector. Instead we calculate directly the 
typical matrix elements of the particle-impurity interac- 
tion within the heavy hole subspace with the spin-orbit 
renormalized wavefunctions. We leave the details for Ap- 
pendix |5] and state the results here: the diagonal terms 
are 

|(tf±3/ 2 |ffF|*± 3 / 2 )l 2 * (3/4)7(1 + l)((3/v ) 2 /N, (29) 

where we neglected small contributions of the spin-orbit 
coupling. An important difference to an analogous result 
for the electrons, Eq. (|28j) is the energy scale. Here the 
typical value for the diagonal Overhauser field is several 
meV, which corresponds to huge external fields of many 
Tesla. The energy splitting of the hole is thus dominated 
by the Overhauser, rather than Zeeman, field. On the 
other hand, the off-diagonal element is non-zero only in 
the presence of the spin-orbit coupling, 

|(*- 3 / 2 |/7f|*3/2)| 2 ~ 2I(I + l)(X (3/vo) 2 /N. (30) 

The impurity spins may induce transition (precession) of 
the heavy hole spin due to the transversal component of 
the Overhauser field, which is smaller by a factor of Aq 



compared to the diagonal component. For our param- 
eters the transversal component is of the order of tens 
of /xeV, so for the heavy hole spin precession to occur, 
the two spin opposite states have to be degenerate with 
respect to this energy [which normally does not occur, 
because of the diagonal term, Eq. ([28])]. 

Having compared the typical energy splittings of the 
particle induced by the effective Overhauser field, versus 
the external magnetic field, we are now ready to discuss 
the crucial assumption for the derivation which will fol- 
low. It is the assumption that the particle is fixed to 
its ground state by an energy gap, irrespective of the 
evolution of the impurity ensemble. This requires that 
spin flips of impurities cost much less in energy than the 
particle transitions 

AEin < AE p . (31) 

For electrons, this assumption is guaranteed as both the 
particle and impurities spin flip costs are dominated by 
the Zeeman energy, proportional to the magnetic mo- 
ment, which is much larger for the electron than for a 
nuclear spin, /ii mp = /ijy ~ W 3 ^b- On the other hand, 
for holes for which the particle and impurity magnetic 
moments are comparable, the above condition is also ful- 
filled since the particle spin flip energy cost is dominated 
by the Overhauser field. 



E. The effective Hamiltonian 

Once the particle is fixed to its ground state (the sub- 
stantial gap assumption), the particle excited states can 
be integrated out perturbativel y 50 ' 66 ' 67 resulting in an ef- 
fective Hamiltonian for the impurity ensemble H c g. For 
this purpose we split the interaction Hamiltonian to 

Hp = Hp + H' F , (32) 

where the diagonal part, 

H% = (* G |77 F |* G )=^K»-P, (33) 

n 

together with the external field, defines the unperturbed 
Hamiltonian Hq — H p + H n z + H F and the basis, so that 

(^ p (g>l A \H nZ + H F \^ q (g)2 B } oc 5 pq 5 AB , (34) 

where T A , 1b denote arbitrary basis states of the impu- 
rity ensemble. We also note that 

(* G ®1 A \H' F \^ G ® X B ) =0. (35) 

Using Lowdin theor y 68 ' 69 the matrix elements of the 
effective Hamiltonian, in the lowest order in the non- 
diagonal part, H' F , are 
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(l A \H cS \l B ) = {^ G ®1 A \ H a 



£ 



1/2 



E, 



GX A 



E, 



1/2 



E GX B 

I — 



E 



pi' 



)T)(y p ®r\H' F \* g ®Ib), (36) 



where the summation proceeds through the excited par- 
ticle states and a complete basis of impurities. We now 
use the substantial gap assumption, which states that all 
system states reachable by the interaction H' F have the 
energy dominated by the particle energy, so that we can 
approximate Eqi — E p z* ~ Eq ~ E p . By this relation 
the summation over the impurities gives an identity, 



{Ta\H cS \T b ) 



>T A \H n 



\%){*p\ 



p^G 



Er 



E r 



H' F \^ G ®1 B ). 



(37) 



Since the impurity states now only sandwich both sides 
of the equation, we can equate the operators 



H eS = (^ G \Ho\^ G ) 



p^G 



(* G \H' F \* p )(* p \H' F \v G ) 

Eg - E p 



(38) 

Even though this looks like the standard second order 
perturbation theory result, note that even after taking 
matrix elements with respect to the particle states, the 
expressions still contain the quantum mechanical opera- 
tors of the impurity spins. On the other hand, by taking 
the expectation value, the particle degrees of freedom dis- 
appear from the effective Hamiltonian. The first term is 
a sum of the particle ground state energy and the impu- 
rities energy in the Knight field, 



(* G |ff |*G> =E G + J28 n -I n - 



(39) 



The first term defines the spin flip energy cost and the 
spin quantization axis, according to Eq. (|21[) . The in- 
teractions described by the second term can be classified 
as spin-conserving (spin-non-conserving) according to ro- 
tated operators I components parallel (perpendicular) to 
a global axis §o, as we will show below. To further demon- 
strate the usefulness and generality of Eq. (|43|) , we show 
that known results follow as special limits, and how the 
consequences of the spin-orbit coupling on the impuri- 
ties interactions can be drawn from the formula. We also 
note that the derivation would proceed in the same way 
even if G were not the particle ground state. The only 
requirement for the validity of Eq. (|4"3")l is that the state 
G is far enough in energy from other particle states so 
that Eq. (|31l) is valid. For example, thermal excitations 
of the particle would result in a thermal average of the 
effective Hamiltonian (the vectors A and energies B do 
depend on G). We do not pursue a finite temperature 
regime further here, and assume the thermal energy k F T 
is small such that the particle stays in the ground state. 

Before we evaluate vectors A in specific cases, we 
note an important property of the effective Hamiltonian. 
Namely, for both holes and electrons, the lowest excited 
state is much closer to the ground state (split by the Zee- 
man energy) compared to higher excited states (split by 
orbital excitation energies) . If the mediated interactions 
are dominated by this low lying excited state, we can 
write 



To simplify the notation of the second term in Eq. (l38l) , 
we introduce 



(*G\Hp\* p ) = (* G \H F \%) = ]TA" F 



(40) 



so that the p-state dependent complex vector A is 

A B = - J 9<* G |<5(R-R n )J|* p ). (41) 

We now transform vectors A and spin operators I into the 
coordinate system along the total field of each impurity 



A n = j£t[A], r = R B l[i] 



(42) 



The z-component of a rotated vector is its projection 
along the direction of the local total field, e.g. I z = I ■ z' 
and similarly for A. Omitting the constant Eg, we 
rewrite Eq. (|38[) with the new notation and arrive at our 
main result 



i 



p^G n,m 



E G -E, 



■(A"-P)(A m -I m ) t . 

(43) 



H cS = 5> n/ "+£ ^^(A^P)(Aft-r»)t, (44) 



where we denoted G =f, the sum was restricted to a sin- 
gle term p =],, and A"^ is the vector corresponding to 
the ground/excited state being the spin up/down state. 
From the relation A™^ = (A"^)t it is straightforward to 
see that had we started with the particle being in the 
spin down state G =1 and restricted to the closest state 
in the spectrum p =f, the mediated interaction Hamilto- 
nian would be the same as in Eq. (I44[) , but for a minus 
sign in the second term, coming from the denominator. 
This crucial property, which results in the particle spin 
decoherence being to a large extent removable by the spin 
echo protocols ) 15 ' 16 is thus preserved in the presence of 
the spin-orbit coupling. 
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Effective Hamiltonian symmetry and magnitude 
of the spin-non-conserving interactions 



(45) 



For electrons, we get from Eq. ([TH)) 
A" = R B l [A"] = e£ R B l o R Un [(J G \ J \J P 



where we denoted the position dependent energy 

£p = -P<S>* G (Rn)%(Rn) ~ -P/V. (46) 

Consider first that the magnetic field is small such that 
the total effective field in Eq. (|20|) is dominated by the 
Knight field, rather than the external field. In other 
words, the local impurity quantization axis is collinear 
with the local particle spin direction. Then = Ru n 
and we get from Eqs. (|43|) and (|45|) the effective Hamil- 
tonian in the following form 



G p jn jrn 

TP Tp Z Z 



(47) 



We have split the summation over the particle excited 
states into those with the same, and the opposite spin as 
is the spin of the ground state, corresponding to the sec- 
ond, and the third term in Eq. (|47p . respectively. Equa- 
tion (|4T|) makes it clear that there is a conserved quantity 
even in the presence of spin-orbit coupling, though it is 
neither the energy nor the total spin along any axis; it 
is the number of impurity spins locally aligned with the 
particle spin, equal to Y] I z . This result is very general, 
as it is based only on the form of the spin-orbit coupling, 
which gives a single unitary operator U for the whole par- 
ticle spectrum. Restricting to the lowest excited state, as 
in Eq. flU), we g et the standard result£££! 



jmjr 



A (48) 



generalized to include the effects of the spin-orbit cou- 
pling. 

For the electronic case, we are, however, more inter- 
ested in a different regime, where a finite magnetic field 
breaks the above discussed symmetry and sets a global 
quantization axis for impurities, so that the Zeeman en- 
ergy dominates in the total effective field in Eq. (|20|) . We 
then have R&™ ~ 1 and Equation (|45|) can be then 

evaluated explicitly, using Eq. 0. Instead, we estimate 
the effects of weak spin-orbit coupling, which guarantees 
that R m -C I so, by expanding the rotation operator up 
to the first order as 



l + 0{r m /l so ). 



(49) 



The pairwise interaction in the effective Hamiltonian 
then appear in combinations such as [see Eq. (|D1|) in 
Appendix ID] 



where 7,7' = O[l/l so ]. This is the most important 
message for the electron case, that the spin-orbit cou- 
pling results in the spin-non-conserving interactions in 
the impurity ensemble, which are, compared to the spin- 
conserving ones, suppressed by a position dependent fac- 
tor of the order of the ratio of the confinement and spin- 
orbit lengths. 

We now turn attention to a hole dot, taking the low- 
est state in the heavy hole subband as the ground state 
G = 3/2,00,0. The closest excited state, which gave by 
far the dominant contribution in the electronic case, is 
the spin opposite heavy hole state p = —3/2,00,0. The 
corresponding vectors A scale as (see Appendix [Dlfor full 
expressions) 



A* 



A'_ 



e p AQ. 



(51) 



To quantify the prefactor in the second order term, v4™ , 
we would have to go to the next order in the perturba- 
tion expansion of the wavefunctions. However, this is 
not necessary as this term does not enter anywhere in 
the subsequent discussion. We conclude from Eq. ([ST]) 
that the spin-conserving interactions mediated by the 
lowest heavy hole excited state are proportional to the 
second power of parameters A [through terms such as 
A-A*_I + IJ\, the same as the spin- non-conserving ones 
[e.g. A_A* Z I + I Z ]. This is a drastic difference to the elec- 
tron case, where the spin-conserving interactions domi- 
nate. 

Let us now consider the light hole subband. Taking 
p = 1/2, 00, 0, we get (see Appendix |Pl) 



Al 



A n 



e p A 0, A z ~ e p Ai 



(52) 



(50) 



The light hole excited state does mediate spin-conserving 
impurity interactions [through A+A+I—I+]. Com- 
pared to these, the leading spin-non-conserving term 
[A + A* Z I-I Z ] is suppressed linearly in A. The energy de- 
nominator in the effective Hamiltonian is of the order 
of 100 meV for the light hole states (typical light-heavy 
hole band offset) versus a few meV offset of the lowest 
heavy hole excited state. For our parameters, this energy 
penalty is almost exactly compensated by much larger 
matrix elements for the spin-non-conserving interactions 
and more than compensated for the spin-conserving ones. 
We conclude that the spin-alike light hole state is the 
most efficient mediator of the spin-conserving interac- 
tions in the impurity ensemble, and rather efficient medi- 
ator of the spin-non-conserving ones. As a direct conse- 
quence, and unlike for electrons, the decoherence induced 
by the hole mediated evolution of the impurity bath will 
not be removed by the hole spin echo. 



IV. PHONON INDUCED SPIN RELAXATION 
OF THE IMPURITY BATH 

We now use the results of the previous section to cal- 
culate how fast the impurity ensemble spin relaxes. The 
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first and the second term of the effective Hamiltonian, 
Eq. (|43p . induces flip of a single impurity and a pair of 
impurities, respectively. For the former, terms with in- 
plane components of I, while for the latter terms such 
as the last two terms in Eq. (l50l) are required for spin- 
non-conserving transitions. As the initial and final state 
energies differ, in general, we consider that the transi- 
tion is assisted by phonons, which provide for the energy 
conservation. 

We consider several possible mechanisms how phonons 
can couple to the impurity bath and make order of mag- 
nitude estimates for the resulting relaxation rates. We 
find that the most efficient relaxation is due to the piezo- 
electric field spatially shifting the particle, leading to a /is 
relaxation time for Mn spins. It is known that phonons 
are ineffective in relaxing nuclear spins j2i still we eval- 
uate the resulting rates also for electrons, as we treated 
electrons and holes on the same footing, the derived for- 
mulas apply for both. We find a 10 11 s relaxation time 
for nuclear spins. 



A. Particle-phonon interactions 

The phonon-impurity interaction Hamiltonian Hi is in 
general a function of the local lattice deformation due to 
the presence of acoustic phonons, 

Q A V 

Here the phonon wavevector is Q, polarization is A (one 
longitudinal A = I and two transversal ones A — tip) 
with eQ>, a real unit vector (eQA = — e_Q,\), Vo is the 
crystal volume, p is the material density, Hujqx = fic\Q 

is the phonon energy, c\ is the phonon velocity, and <Zq A 
is the phonon creation operator. 

In a polar material, such as GaAs, the lattice deforma- 
tion is accompanied by a piezoelectric field, which is the 
gradient of the following potential 

X (Q x Qy{eQ\)z + QzQx{eQx)y + QyQ z( e Q\)x"j , 

(54) 

with S the piezoelectric constant. 

In addition to the previous, the deformation of the 
lattice shifts the electronic bands, quantified as the de- 
formation potential Vdp = —<J div <5R, which thus reads 

^ = 'E Jj^ Q ciQ R fa + a -Q<) ' ^ 

Q ' 

with a the deformation potential constant. 



In what follows we will see that a relative shift of the 
impurity and the particle, which we denote by d, will 
induce impurity-phonon coupling, leading to the impu- 
rity spin relaxation. Since impurities are tied to atoms, 
the phonon displacement is obviously such a relative 
shift d = <5R, which we call "geometric". However, 
the phonon induced electric fields E also lead to shifts. 
Namely, adding the potential of an in-plane field to that 
in Eq. ([2]) amounts to a shift of the quantum dot position 
by d = e~El 2 /tiui (electrically induced shifts along the per- 
pendicular direction are much smaller, as the wave func- 
tion is much stiffer along z due to stronger confinement). 
If the particle follows these potential changes adiabati- 
cally, which we assume, such a shift is equivalent to the 
shift of the impurities, which are fixed to the lattice, by 
— d. Since the phonon electric fields are proportional to 
the displacement 5R, we can write a general expression 

|d|~a|5R|, (56) 

with a dimensionless factor a. For the geometric shift 
mechanism a = 1 by definition. For the piezoelectric 
field, comparing Eqs. (|53")l and (|5H) . we get 

a = 2Ql^. (57) 

Finally, the deformation potential gives 

a = a{Ql) 2 /huj. (58) 

We specify the dimensionless factor a in Table fl] As it 
enters the relaxation rates in the second power (as we 
will see), we can immediately quantify the relative impor- 
tance of the three considered channels. Piezoelectric field 
is the most effective, for both electron and hole cases, in- 
ducing shifts almost two orders of magnitude larger than 
the geometric shift. The electric field from the defor- 
mation potential is comparable to the piezoelectric for 
holes, and much smaller for electrons, which are deeply 
in the long phonon wavelength limit, Ql <C 1. We note 
that the geometric shifts will be in fact somewhat more 
effective than it seems from the table, as they may (un- 
like the electric fields) shift the wavefunction along the 
perpendicular direction. This results in an effective en- 
hancement of a by a factor of irl/w, which, however, is 
not large enough to change the order of importance fol- 
lowing from the Table. 

We will calculate the relaxation rate T due to a general 
phonon-impurity interaction Hi by the Fermi's golden 
rule. For a given phonon polarization A it reads 

r = Y E \(I'mi)\ 2 S(Ex - E x , - hu QX )N Q , (59) 
Q 

where X and I' denote the initial and final state of the im- 
purities, and we are interested in transitions where these 
two states differ in spin. The phonon occupation factor 
Nq = uq + 1, and Nq — hq, if the energy of the initial 
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a 


piezoelectric deformation geometric 


electron 
hole 


46 0.0038 1 
38 17 1 



TABLE I: Values for the dimensionless coefficient a, the ratio 
of the induced impurity shift versus the phonon displacement 
for various shift mechanisms (columns) and particles (rows). 
Parameters from Appendix [E] were used (GaAs and ZnTe for 
the electron and hole case, respectively), the phonon wavevec- 
tor for electronic case was specified choosing B = 1 T. 



state is larger, and smaller than the final state, respec- 
tively, with riQ = l/[exp(hujQ\/kBT) — 1]. The energy 
conservation fixes the phonon wavevector magnitude to 
\E X — E x <\ = hc\Q, by which we get 



r = -^iv Q |(i'|i^>l 2 - 



(60) 



Trh 2 c\ 

The overline denotes the angular average, 

7(Q) = (l/47r) fdSlf(Q), (61) 
over directions of the phonon wavevector Q. 

B. Spin-phonon coupling due to impurity shift 

Assuming the shifts are small, we get the impurity- 
phonon coupling as 



E dH c ff 
d--^HH-H. 



<9R 



(62) 



To calculate the spatial derivative of the effective Hamil- 
tonian, Eq. (l43l) , it is easier to first evaluate the deriva- 
tive of vectors 3 and A in the original coordinate system 
which does not depend on the position and then to trans- 
form them into the locally rotated coordinates. 

Let us start with the first term of the effective Hamil- 
tonian, Eq. (|43|) . The finite derivative of the total field, 
Eq. (|2"0"]l. is due to the spatial dependence of the Knight 
field, 



(d • dn)B n = (d • 9 R )K" 



(63) 



transversal components of which give the spin increasing 
transition rate for impurity n as 



r (i) = 



VoQ 2 

irh 2 c\ 



N Q [r'_ ■ (d ■ a R )K]2|(r +1 |^|/")| 2 , (64) 



and an analogous expression follows swapping the sub- 
scripts ± for a spin decreasing transition. Our goal is 
an order of magnitude estimate for the rates as the one 
in Eq. ([64]). which allows us to simplify: We replace the 
spin dependent matrix elements of the raising/lowering 
operators by 



" ' ^ ' I U 1 1 71 ) | ^ 



/(/+l)-/"(/"±l)~/ 2 , (65) 



choose the direction of the phonon wavevector that gives 
the highest contribution, instead of performing the an- 
gular average in Eq. (|6Tj) . and denote 



VK_ = <TV_ 



Finally, we use Eq. (|56|) and SR 



■ 0r)K. (66) 
y / h/2Vopc\Q to write 



r (i) 



a 2 I 2 N Q Q 
2ttHc 2 x p 



(VK_ 



(67) 



a general form for the relaxation rate estimate, which we 
a few lines below evaluate for specific cases. 

Let us start with the electronic case. The transition 
energy is dominated by the external field hc\Q ~ \gpNB\ 
and is much smaller than the thermal energy, so that 
Nq ~ k B T '/hc\Q. We evaluate the derivative of the 
Knight field in Appendix ICl see Eq. (|C6|) . getting 



VK± ~ {P/V)l-\ 
For the dominant piezoelectric mechanism we get 
9I 2 k B T(gp N B) 2 l A 



r (i) 



2tt 3 



hE 2 pz 



w 2 l 



2/2 



(68) 



(69) 



where the energy E pz — y/ffc^pfEmfi is a material con- 
stant. Evaluating parameters of GaAs, we get a mi- 
nuscule rate T^ 1 ) ~ 2 x 10 _11 s _1 , choosing transversal 
phonons, external field 1 Tesla and temperature 1 Kelvin. 

We now turn to holes. The transferred energy is now 
given by the Knight field, rather than the external field, 
\Ei — Ex* | ~ J(3/V and we again consider a high temper- 
ature limit, UbT > \Ex — E%' |. As we show in Appendix 
[Cl the formula in Eq. is changed into 



VK± ~ tf/V)(V3\i/l), 



(70) 



showing that quantities of the form l/X can be seen as 
an effective spin-orbit length for holes. Equation (|69|) 
can be then used putting for the "spin-orbit length" the 
one just described and replacing the external Zeeman 
energy by the Knight field. For convenience, we give the 
relaxation rate explicitly 



-(i) 



3 5 / 2 J 2 k B T (3 2 X 2 



2 5 7T 4 



H E 2 w 4 l 2 ' 

pz 



(71) 



which for ZnTe parameters and temperature 10 K yields 
a modest rate 0.35 /is -1 . The rate is second order in 
the "spin-orbit strength" and, unlike for electrons, grows 
very fast upon making the dot smaller, since now the 
effect of Knight field, inversely proportional to the dot 
volume, dominates over the effect of the shift being larger 
for softer dot potential. 

Comparing the two terms of the effective Hamiltonian, 
Eq. (|43t , and using the results of Appendix [Dl the pair- 
wise spin transition mediated by the lowest heavy hole 
excited state relates to the single spin-flip rate by 



at (2) , /r« 



V/ 2 f 2- A ° 



p/v 



Ai Ei — E* 



(72) 
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which can be evaluated as 10 2 . Similarly, we get for the 
mediation through the light hole state 



JYp(2) /p(l) 



N 



4Ax A lh J 



(73) 



where a much larger energy offset, A;^, is partially com- 
pensated by a larger matrix element. Note that the num- 
ber of pairs available for a flip is of the order of N times 
larger than the spins themselves, so that when compar- 
ing the first and the second order rate, the latter should 
be multiplied by N, what we did in the previous two 
equations. 

Just for completeness we note that for electrons a sim- 
ilar relation between the first and second order rates 
holds, 



-(2) 



E± - R 



r <r 



(i) 



(74) 



but the ratio is much smaller, at the external field of 1 
Tesla by 10 orders of magnitude. 

We now consider additional mechanisms of the 
phonon-impurity couplings, through which impurity spin 
relaxation may arise. 



C. Valence band shifts 

The phonon induced lattice compression changes the 
bandstructure - the bands are shifted. Shifts different for 
the bands of the particle ground and mediating excited 
state p result in the impurity-phonon coupling through 
the second term of Eq. (|43|) , by changing the denominator 
by AVdp- As the two states have to belong to different 
bands, such a coupling may arise only for the case of 
holes and takes the form 



Hi 



AV, 



DP 



H 



effj 



Hh 



(75) 



expanding the effective Hamiltonian up to the lowest 
order in the band shift difference AVdp = — (chh — 
aih)dxvSH, which is the difference of the deformation po- 
tentials for the heavy and light hole valence bands. We 
can thus relate the band shift mechanism to the position 
shift one, comparing Eq. (|62p with Eq. (|75[) . We find that 
the latter is described by an effective constant a 



ote$ = (om - vih)crQl/A 



Hi- 



(76) 



If we estimate the potential difference by the typical value 
of the potential itself, (<Jhh — &lh) ~ 5 eV, which is likely 
an overestimated value, we get a e g = 17, so that the cou- 
pling through band shifts leads to a rate at most com- 
parable to (and most probably much smaller than) that 
described by Eq. ([75]). 



D. Renormalization of the spin-orbit length 

Phonon induced renormalization of band offsets influ- 
ences the spin-orbit couplings. This is evident from the 
expressions for the coefficients A which are inversely pro- 
portional to the light- heavy hole offset Aih, see Eqs. (|A9I) . 
Even though the coupling in the form of the constant 
a is described by the same formula as in Eq. ([75]) . the 
substantial difference is that now the phonon-impurity 
coupling arises also through the first term of the effective 
Hamiltonian because fluctuations in spin-orbit fields in- 
duce fluctuations in the Knight field. The corresponding 
a e g is the one given in Eq. (|76l) and the relaxation rate 
is that in Eq. (|67|) . so that it does not exceed the rate 
due to the piezoelectric shift mechanism. For the case of 
electrons, the effective constant follows in an analogous 
form, 



a = (<x e - a h )aQl so /A, 



(77) 



where A is of the order of the conduction- valence band 
offset, which enters the definition of the spin-orbit cou- 
plings l/l so . The numerical value for a is much less than 
one even taking (cr e — ah) ~ cr e , so that this channel is 
negligible with respect to the geometrical shifts. 



E. Phonon induced spin-orbit interactions 

Finally, we estimate the influence of spin-phonon cou- 
pling through an additional spin-orbit interaction aris- 
ing in the presence of a phonon-originated electric field. 
We assume the new spin-orbit interaction strength re- 
lates to the one we considered in Sec. II, referred to as 
"old" , through the ratio of the internal (interface) elec- 
tric field Ei nt and the phonon induced electric field, the 
latter given as the gradient of the appropriate poten- 
tial, Eqs. (HU) or (f5"5l) . Since the phonon- induced spin- 
orbit interaction ("new") arises from an electric field, it 
is of the Rashba functional form, and one can take its 
effects to be additive to the "old" one. For electrons, 
this means U tot = U i d U new « (I + o i d )(I + O new ) « 
I + (Odd + O new ) , which amounts to additive inverse of 
the spin-orbit lengths l// so — > l/l so + l/l P h-so- By in- 
specting Eq. (|C5I) . we estimate the effective interaction 
to be described again by Eq. (|67l) with the constant 



BQl 

e-Eint 



(78) 



for the piezoelectric phonon field (a much smaller defor- 
mation field corresponds to the numerator replaced by 
crQ 2 l). The numerator evaluates to 10 6 V/m, which is 
not supposed to be much higher than the internal field, 
so that again, we find that this mechanism is less im- 
portant than the one due to the piezoelectric shift. We 
come to a similar conclusion for holes (though we do not 
show calculations in detail here). Namely, even though 
the phonon induced piezoelectric field is much stronger, 
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of the order of 10 8 V/m, the spin-orbit terms in Kohn- 
Luttinger Hamiltonian are less effective in inducing light- 
heavy hole mixing than the terms we considered explic- 
itly in previous sections (see Appendix [F]). 



V. CONCLUSIONS 

We have analyzed the interactions within an ensemble 
of impurity spins, which are mediated by a confined spin- 
orbit coupled particle. We have considered two physical 
systems where such mediated interactions are of great 
importance: III-V (GaAs) electronic lateral quantum dot 
with the impurities being nuclear spins, and II- VI (ZnTe) 
self- assembled Mn-doped quantum dot populated by a 
hole. Our focus has been on the consequences of the 
spin-orbit coupling on the character of the mediated in- 
teractions. 

We have derived an effective Hamiltonian for the impu- 
rity ensemble, treating the particle-impurity interaction 
perturbatively. The form of this Hamiltonian allowed us 
to quantify the degree to which the conservation of im- 
purities spin is broken in the presence of the spin-orbit 
coupling of the particle. We have found that for the elec- 
tron case, the spin-non-conserving terms are suppressed 
relative to the spin-conserving ones by a small factor, the 
ratio of the confinement length and the spin-orbit length. 
The lowest electron excited state is the most effective me- 
diator, what results in a decoherence being removable by 
the electron spin echo even in the presence of the spin- 
orbit coupling. In the case of holes, the spin-conserving 
interactions are most efficiently mediated by the lowest 
light hole state, while the spin-non-conserving ones by 
the lowest heavy hole state. The induced decoherence is 
then not removable by a hole spin echo anymore. 

As a direct application of the derived effective Hamil- 
tonian, we have calculated the rates of a phonon-assisted 
impurity spin relaxation, which arises only in the pres- 
ence of the spin-orbit coupling in the particle Hamilto- 
nian. We have considered several coupling mechanisms, 
by which impurity spins couple to phonons. We have 
found that the most effective is the piezoelectric field in- 
duced shift of the particle wavefunction, with a typical 
relaxation time of 1 fj,s in a 10 nm self-assembled strain- 
free quantum dot. The rate grows upon making the dot 
smaller, or diminishing the heavy-light hole splitting by 
strain, possibly to nano seconds for reasonable dot pa- 
rameters. 

While we have focused on a single particle, our anal- 
ysis of the spin-non-conserving mechanisms already pro- 
vides insights in the magnetic ordering in quantum dots 
with multiple occupancy. For example, the spin-non- 
conserving mechanisms are the key in understanding the 
prediction of piezomagnetic quantum dots or, equiva- 
lently, nonlinear magneto-electric effects.- 7 "^ Changes in 
the shape of lateral confinement (from circular to ellip- 
tical) controlled by the pair of the gate electrodes have 
been demonstrated experimentally to alter the particle 



configuration from vanishing to finite spin in nonmag- 
netic quantum dots^ This principle provides intrigu- 
ing possibilities for the control of magnetic ordering in 
dots with added Mn impurities where such changes in 
the particle spin, through exchange interaction, would 
reversibly control the magnetic ordering of the nearby 
Mn spins Q The characteristic time scale for the related 
magnetic polaron formation, similar to the better studied 
quantum well structures, should be largely determined 
by the anisotropic spin-spin interactions of impurities 
which explicitly do not conserve the total spin of Mn 
atoms^SiZi 

Beyond epitaxial dots that we have presently exam- 
ined, recent experimental advances in colloidal quan- 
tum dots warrant also future considerations. Typi- 
cally they are easily synthesized II- VI materials, such 
as ZnTe, ZnSe, CdS, and CdSe^^ which offer a large 
size- induced tunability of the transition energies and long 
spin decoherence times XL Magnetic doping 78 of these col- 
loidal dots provide an opportunity for a versatile control 
of magnetic order as well as lead to robust magnetic po- 
laron formation with effective internal magnetic field ap- 
proaching 100 Ti^rH 
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Appendix A: Kohn-Luttinger Hamiltonian 
perturbative eigenstates 

Here we derive hole eigenstates in the lowest order per- 
turbation theory. We neglect the influence of the conduc- 
tion and spin-orbit split-off subbands and consider only 
the light (J = ±1/2) and heavy (J = ±3/2) holes in 
the Kohn-Luttinger Hamiltonian^ Hjji . The diagonal 
elements are 

h 2 

#±3/2,±3/2 = [£'(71 - 2 72 ) + {K + fc')( 7l + 72 )] , 

h 2 

H ±1 /2,±x/2 = [Khi + 272) + {K + ^)( 7 i - 72)] , 

(Al) 

with Kk = P the hole momentum operator, mo the free 
electron mass, and 71, 72, and 73 (below) the Luttingcr 
parameters. Together with the in-plane V(r) and het- 
erostructure V z (z) confinement potentials, assumed to be 
those in Eq. ^ , the kinetic terms in Eq. (jAll) define the 
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dot unperturbed eigenstates (normalization omitted) 

< m ,k = r^e-^^L^iryi^e^siniknz/w). (A2) 

We used cylindrical coordinates (r, <j), z), L™ are the asso- 
ciated Laguerre polynomials, lj the in-plane confinement 
length (which differs for heavy and light holes due to their 
different masses) and w the heterostructure width. The 
Fock-Darwin states are labelled by the principal and or- 
bital quantum numbers n and m, and k labels excitations 
in the perpendicular potential. The corresponding ener- 
gies are 



Ej. 



'nrn.k 



HQj(2n 



1) 



r^ 2 , (A3) 



where the in-plane excitation energy is parameterized by 
the mass and confinement length, 



mjlj 



(A4) 



The in-plane masses are given by Eq. (|A1[) : m-1-3/2 = 
m hh = m /(7i +72), rn ±1/2 = m Ul = m /(7i -72)- We 
parameterize the in-plane electrostatic potential choos- 
ing a certain value for the heavy hole in-plane excitation 
energy £3/2,01,0 " £3/2,00,0 = T&lhh, which then specifies 
the confinement lengths. The light hole excitation energy 
and confinement length 

ftlh = Qhh{mhh/mih) 1/2 , hh = lhh{m hh /mi h ) 1/4 , 

differ from the corresponding heavy hole quantities due 
to a different in-plane mass. The energies of the hard- 
wall eigenstates also differ for heavy and light holes due 
to different out-of-plane masses, which are 7710/(71 + 272), 
and mo/ (71 — 272), respectively. We set w by choosing a 
certain value for the light- heavy hole splitting £3/2.00.0 — 
£1/2,00,0 = A lh . 

Taking the above eigenstates as the basis, we now per- 
turbatively take into account the off-diagonal elements of 
the Hamiltonian, 



£±3/2, ±1/2 — ± 



2t7Jo 

h 2 



2v / 37 3 /c T fc z , 



£±3/2, T i/2 = 7, [j2(kl - k 2 y ) =F 2i lz k x ky\ , 



(A6) 



2mo 

£±3/2,^3/2 = = £±1/2,^1/2- 



We employ the non-degenerate perturbation theory 



|*Ja) = \J) ® |^) 



J' a'^Ja 



(* j a ;\Hj,j\$ j a ) 

E Ja - Ej, a . 



I J') 



(A7) 

where we use the notation introduced below Eq. (QJ, so 
that a includes two quantum numbers of the in-plane 
Fock-Darwin state and one of the perpendicular hard- 
wall state. To proceed, we neglect high energy excitations 
n > and k > 1 and adopt the axial approximation 

h 2 



H 



±3/2,^1/2 



2mo 



V3 7 fc|, 



(A8) 



with 7 = (72 + 73)/2. With these simplifications, the 
otherwise infinite sum for \^ j a ) simplifies to only a single 
term for each J' ^ J and can be given explicitly.—^ We 
finally get Eq. © with the admixtures 



Ai 
A 



7 3 a/3< 



m^lihW A* + Ml 



Hi 



7 



3/2k 



m l 2 A lh + 2Ml Ul ' 



(A9a) 
(A9b) 



where A* = £1/2,00,1 
hole offset and 

«=(<oKo>- 



£3/2,00,0 i g the z-excited light 



(mfch/mih) 1 /* + {mih/mhh) 1 '^ 

(A10) 

is the ground state overlap, which differs from one due to 
the in-plane mass difference. Finally, the dimensionless 
matrix element £ is defined by £ = — ku(l|fe z |0) =8/3. 

Taking a heavy hole 3/2, the admixture of 1/2 light 
hole scales as 1/wl, costs the in-plane plus perpendicular 
orbital energy (the latter is even larger than the light- 
heavy hole offset) and leads to an admixture with a very 
different z-profile (compared to the main wavefunction 
component). The admixture of the -1/2 light hole has a 
smaller numerator, proportional to l/l 2 , but costs only 
the in-plane orbital energy (several times smaller than 
the light-heavy hole offset) and has the same z-profile as 
the main component. 

Along the same lines we get Eq. (fTTj) with 



A'i = 



An 



73 



V3k£ 



molhhw A*' + hflhh ' 

h 2 -i^/W^k 



m l 2 A L 



2Mi 



1 1 h 



(Alia) 
(Allb) 



l hh 

where A*' = £3/2,00.1 — £1/2,00,0 i s the z-excited heavy 
hole offset. 

For completeness, we list the hole Zeeman term§£ 
H hz = 2 K [i B J ■ B + 2 qi i B £ JfBi, (A12) 

i=x,y,z 

which can be written for the heavy hole subspace as 

H hh , z = g hh n B [J/3] z -B z , (A13) 

with J/3 = er/2 the pseudo spin operator and g^h = 
6k + 27<7/2 « 2 for GaAs^ 

Appendix B: Spin matrix elements for holes 

Here we calculate the matrix elements of the spin op- 
erator between perturbative eigenstates of the hole. For 
the purposes of this appendix, we shorten the expres- 
sion in Eq. ([3]) introducing $ = $00,01 a = Ai$ i,i anc ^ 
b = Ao$&,„ to 

1*3 



I*- 



> 3/3/ = *|3/2) + a|l/2) 
. 3/a ) = **| - 3/2) - a*| 



6|-l/2), (Bl) 
- 1/2) +6*|l/2>- (B2) 
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We denote J± = J x ± U y , so that J x = (J + + J-)/2 
and J y = (J + — J_)/2i and since the orbital operator in 
all the matrix elements is the delta function 5(R — R n ), 
all the complex amplitudes below should be evaluated at 
the position of the particular impurity, e.g., $ — > $(R„). 
Listing only the leading order in small quantities Ao,i, to 
which a and b are proportional, we have 

(*3/ 2 |J.|* 3 /2> = (3/2)|$| 2 + 0(A 2 ), (B3a) 
(*3/ak+l*3/ 2 ) = V34>*a + 0(\ 2 ), (B3b) 
(* 3 /2|J-|*3/ 2 ) = V3$a*+0(A 2 ), (B3c) 

from where Eq. (fl9|) follows directly. The time reversal 
symmetry gives (using (Ta\b) = —(a\Tb)*, and T 2 = -1) 

(*-3/a|J|*-3/2) = -(*3/2|J|*3/2), ( B4 ) 

so that the spin expectation value changes sign upon in- 
verting the hole spin. To evaluate the off-diagonal ele- 
ment of the Overhauser field and the vectors A we need 

(*_ 3/2 |J z |* 3/2 ) = ab, (B5a) 

(*-3/ 2 k+|*3/2) = 26 2 , (B5b) 

(*_3/ 3 |J-l*3/a) = 2v / 3$& + 0(A 2 ), (B5c) 

from where Eq. (jBTJ)) follows. 

Analogously we define short hand notations for the 
light hole like wavefunctions as $' = $oo. , a ' = ^'i^oT 1' 
and b' — Aq$^ , to write 

|* 1/a ) = $'|l/2)+a'|3/2)-6'|-3/2), (B6) 
|*-i/a) = *'*|-l/2)-a'*|-3/2)-6'*|3/2).(B7) 

The light hole mediated effective interaction is given by 

(*3/2|J*|*i/ 2 ) = (3/2)$V + (l/2)a*$', (B8a) 
(*3/a|J+l*i/2> = V3$*$' + 0(A 2 ), (B8b) 
(* 3 /2|J-l*i/2> = 2$V+0(A 2 ). (B8c) 

Finally, the coupling through the spin opposite light hole 
state are given by 

(*3/2|J*|*-i/ 2 ) = -mWb'* - (l/2)6*$*'CB9a) 
(* 3 /a|J+l*-i/2) = 2a*$'*+0(A 2 ), (B9b) 
(*3/a|J-|*-i/ 2 ) = -V2>b'*a*. (B9c) 

Appendix C: Knight field and rotated coordinates 

In this appendix, we evaluate the Knight field, its spa- 
tial derivative and the corresponding locally rotated co- 
ordinate frame, for the case of electrons and holes. 



1. Electron case 

The Knight field of an electron in the ground state is 

K = -(/?/2)|<& G (R)| 2 (s(R)>, (CI) 

see Eq. (fT8|) . In deriving that, we used the identity 

UJU ] =R U 1 [3] : (C2) 

where the unitary operator of spinor rotation 

U ee exp(-in- 3(f)), (C3) 

corresponds to a three dimensional rotation around the 
unit vector n by angle 4>, 

Ru =exp(-in-10), (C4) 

with (Ikjmn = — iffcmra- 

The spatial derivative of the Knight field, which gives 
the matrix elements for the impurity spin flips in the 
"first order" rates, follows from Eq. (|C1|) as 

VK= [Vln(|$ G (R)| 2 )]K + (Vn so ) xK. (C5) 

The two terms correspond, respectively, to the position 
change in the magnitude and direction of vector K. 

For electrons we consider a regime in which the total 
field is dominated by the external field. The impurity lo- 
cal coordinate system then coincides with the coordinate 
frame defined by the external magnetic field, « 1 
and %' = sq. This allows us to estimate the components 
of VK transversal to the local coordinate system as 

r'± -VR-cr^so)^, (C6) 

which originate in the first term of Eq. (|C5[) and where 
the length 5 — I for a phonon with in-plane polarization 
vector (e _L z) and 6 = w/tt for a phonon with a po- 
larization vector along the growth direction (e || z). The 
second term of Eq. (1C5|) gives a contribution at most as 
large as the first term, or lower, depending on the phonon 
polarization. 

2. Holes 

The Knight field of a spin 3/2 hole, defined by Eq. (fT7|) . 
K" = -f3(* 3/2 \6(R- Kn) J 1*3/2}, (C7) 
follows from Eqs. (|B3[) as 

[K5,K?, iC] = -/3[v / 3Re(a$*),v / 3Im(a$*),3/2|$| 2 ], 

(C8) 

given in the main text in Eq. (|19l) , here using the notation 
from Appendix [0 a = Ai$o'j )1 (R n ) and $ = $^ (R n ). 
Since for holes we are interested in the zero magnetic 
field case (by which $ is real), the local coordinate frame 



15 



is defined by a unit vector along the Knight field, z' = 
K/K, 

z' = [sin</>Re(a)/|a|,sin</>Im(a)/|a|,cos</>]. (C9) 

Here <j> is the angle between th e original an d the rotated 
z axes, cos0 = z • z' = <I> / ^/3 ( E >2 /4 + \a\ 2 . We choose 
the remaining two axes of the rotated coordinate system 
arbitrarily as 



-Im(a)/|a|,Re(a)/|a|,0], 



(CIO) 



and 



y x z 



[cos 0Re(et)/|a|, cos </>Im(a)/|a|, — sin</>]. 

(Cll) 

Using the projectors into the transversal plane of the lo- 
cal coordinate system, r'±, one can evaluate the ampli- 
tudes of the spin-non-conserving terms. As an auxiliary 
result, we note that for any real vector v we have 



r ± • v = v + - 



a* cos 6 ± 1 



a cos cj) =F 1 



\a\ 2 \a\ 2 r 

(C12) 

For example, close to the dot center, the inequality $ ^ 
|a| gives sin0 « |a|/v / 3/2$, cos</> « 1, \a\ ~ Ai<£>, and 
the Knight field derivative follows as 



r ± 



(C13) 



where, again, the length 8 depends on the direction of 
the phonon polarization vector, 5 — I for eQ _L z, and 
8 = w/tt for eQ || z. 



Appendix D: Interactions in the impurity ensemble: 
vectors A and the effective Hamiltonian terms 



The second order Hamiltonian for a given impurity pair 
n, m equals 1/ (Eq — E p ) times the following expression 



(A n • I n )(A ro • l m ) 

rn jm ( a n a m* 

r_ii m (AiA m * 

jnjrn^ A n A rn* 

I n z I™(A n z A™ 
I\I™(A n _A 

r z i i™{A n z A 
ri™(A%A 

IlI™{A n _A™* 



'Z 



m* 
z 



i t + (A m -I m )(A"-I") t 

+ A™Al*)/A+ 

+ A m A n»y 2 + 

+ A , z ,l A n + *)/2+ 
+ A™A n z *)/2+ 
+ A r z n A n *)/2+ 
+ A m Al*)/4+ 
+ A™A n _*)/A 



(Dl) 



If the vectors are expressed in the local coordinate frame 
[that is, all quantities in Eq. (|D1[) with tildes], the first 
three terms are spin preserving (connect states with the 
sum of spin projections along the local spin quantization 



axes), the next four terms change the sum by one (repre- 
senting a single spin flip), and the last two terms induce 
double flips. The complex conjugates arc defined as 

Am* /„/ a m\* Am* ■ Am* I / \ m\* /tyo\ 

A± = (r ± • A ) = A x ^iA y = r T • (A ) . (D2) 

To find the effective Hamiltonian, it remains to eval- 
uate the vectors A. For electrons, we assume that the 
magnetic field dominates the total field for impurities. 
Equation (gSJ) gives (J = 1/2, J' = -1/2) 

A" = -(/3/2)|$ G (R„)| 2 i? [/ „[r'_], (D3) 

where we remind that the vectors x', y', and §o = z' 
form an orthonormal set, with §o along the external field. 
Expanding the rotation operator in the lowest order in 
the spin-orbit length we finally find 



-(/3/2)|$ G (R„)| 2 , Al,A?~0(r n /l)A1, (D4) 



A'! 



from where Eq. (|50[) of the main text follows. 

For holes, we find the vectors A corresponding to the 
interaction mediated by the -3/2 state from Eqs. (|B5j) 
and Eq. ([Ci~2|) as 

r'+ • A = -/% 2 (a/|a|)(cos</> + 1), 
r'_ • A = -V§j8$&(a/|a|)(cos0 + l), (D5) 
t! ■ A = -V3j3<f>b(a/\a\) cos 0. 

The largest spin-non-conserving term of the correspond- 
ing effective Hamiltonian is 

H nm „ fnfm j_ g ^(^2)2^ (pg) 

The vectors A for the spin alike light hole, using 
Eqs. (HU follow as 



r'+ ■ A = -v / 3 / 9$*$'(a7|a|)(cos0+ l)/2, 
r'_ • A = -f3<f>'b*(cos(f>+ 1), (D7) 
z' ■ A = -(/3/2)(3a'$* + a*$') cos0. 

From the above results, we see that the largest spin- 
non-conserving terms in the effective interaction are 



and 



//',":' - J" I'!' [ -^\/3(/3$ 2 ) 2 A 



H nm „ jnjm [ _L ( 3 /4) V^* 2 )^ 



(D8) 



(D9) 



Similarly as before, a differentiation with respect to the 
position, which enters the impurity-phonon rates, brings 
in an additional factor 1/5 in Eqs. (|D6|) - (|D9|) . 
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quantity 


electron hole 


V 
N 

p/v 


3 x 10 4 nm 3 238 nm 3 
1.3 x 10 6 iMnX 5278 
0.13 neV 64 peV 



TABLE II: Effective volume, number of impurities and the 
coupling energy scale. 



Appendix E: Materials parameters 

For the electronic case we assume a GaAs/AlGaAs het- 
erostructure with the following parameters: electron ef- 
fective mass m — 0.067 mo, in-plane confinement length 
/ = 30 nm, quantum well width w = 8 nm, spin-orbit 
length l so ~ 1 pm, electron-nuclear coupling /3 = 4 peV 
nm 3 , material density p = 5300 kg/m 3 , phonon veloci- 
ties ci = 5290 m/s and Cf = 2480 m/s, conduction band 
piezoelectric S = 1.4 x 10 9 eV/m and deformation a = 10 
eV potentials. The g-factors and corresponding energy 
scales are given in Table IIIII 

For the hole case, we list the Luttinger parameters 
71/72/73 of GaAs:SL2£ 7.1/2/2.9, CdTe^.l/l. 1/1.6, 
and ZnTe>22 3.8/0.7/1.3. We take ZnTe as the mate- 
rial of our choice, with p = 5650 kg/m 3 , q = 3550 m/s, 
c t = 2358 m/s, a = 5 eV, 5 = 3.4 x 10 8 eV/m. We 
set the heavy hole orbital energy Mlhh 1° 20 meV, which 
gives Ith — 4.19 nm, lih = 3.82 nm and Mlih = 17 meV. 
We set the light-heavy hole splitting Aih to 100 meV, 
which gives w — 3.24 nm. The hole-impurity interaction 
strength is /3 = 1/3 eV do/4, with ao = 0.61 nm the 
lattice constant. We assume the Mn impurities concen- 
tration is given as iMm the ratio of cation replaced by 
Mn atoms, typically of the order of 1%. 

Properly normalized, the ground state follows from 
Eq. (|A7| as 



<oM,*) = /I sin O 7=7; ex P (-571) . 



from where the quantum dot volume estimate 



(El) 



quantity 


nucleus 


Mn 


electron 


hole 


9 


1.2 


2 


-0.44 


-2/3 


A 4 


fJ,N 


/IB 


\IB 


PB 


I or J 


3/2 


5/2 


1/2 


3/2 


g^i 


57 neV/T 290 /xeV/T 12.7 peV /T 58 ^eV/T 


gpI/k B 


660 pK/T 3.4 K/T 


146 mK/T 


670 mK/T 


B e f[ 


66 peV 


96 peV 


290 neV 


13.8 Xy^ meV 




1.2 mT 


331 mT 


23 mT 


237 zl 2 T 



TABLE III: g-factor, magnetic moment, spin, the correspond- 
ing energy scale and the effective field B c g . For GaAs nuclear 
spins, the g-factor is the average over naturally abundant iso- 
topes. The effective field for impurities is the Knight field, 
B c ff = J(/3/V), for the particle it is the Overhauser field, 
B cB = y/NI(I + l)(p/V). 



the corresponding coefficients A. We find that these al- 
ternatives lead to negligible amount of admixture, com- 
pared to the terms of the Kohn-Luttingcr Hamiltonian 
we considered in the main text. 

In the presence of electric field E there appear the fol- 
lowing term in the hole Hamiltonian (in the notation of 
Ref.ia) 



(Fl) 



Assuming, for convenience, that the field is perpendic- 
ular to the heterostructure interface E = zE this term 
translates into our notation of Sec. [Alas 



H 



±3/2, ±1/2 



±rlf v EV3k T /2, 



(F2) 



where the material parameter rfi we estimate by its 
ZnSe value of -4.1 eA 2 . Taking Eq. (fA7|) for the heavy 
hole ground state ^3/2), the Rashba spin-orbit inter- 



action leads 
|l/2>® 



<f> lh 



^3/2 

to an admixture of the 
> with a coefficient 



V37 



E 



2 A 



ih 



light hole state 



(F3) 



V = 1/ / d 3 R|$(R)| 4 = (4tt/3)w;Z 



(E2) 



gives values in Table [TTJ 

From the above parameters, the hole admixture coef- 
ficients follow as Ao ~ 0.053, Ai ps 0.050, A ~ 0.11, and 
A; ps 0.15. 



Appendix F: Alternative hole mixing mechanisms 

In this section we estimate two additional light-heavy 
hole admixture sources, namely the spin-orbit coupling to 
an external electric field (the Rashba spin-orbit interac- 
tion) and Overhauser impurity field itself. We quantify 
the resulting light-heavy hole admixture by calculating 



To quantify its value, it remains put in the value for the 
electric field E. We estimate it to be of order 10 7 V/m, 
using a very crude electrostatic model; namely we take 
the hole to be a point charge in the center of a half sphere 
of a radius I ~ 5 nm and the electron to be a uniform 
classical charge density a on the half-sphere, which corre- 
sponds to the type Il-quantum dot populated by a single 
exciton. The resulting internal field E = a/2e is smaller 
than the piezoelectric field accompanying the phonon, 
Eq. (|78p . which we estimated to be of order 10 8 V/m. 
Inserting this latter value for E, we find A so ps ^/3 x 10~ 3 , 
so that both the internal electric field as well as phonon 
induced fields induce light-heavy hole admixtures much 
smaller than those we considered in the main text. 

One may also consider the collective field of the impu- 
rities (the Overhauser field) as a source of the light-heavy 
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hole mixing. Using Eqs. (|B8[) we estimate the arising ad- 
mixture of the light hole state 1 1/2} ® |<& q ) into the 
heavy hole ground state 

\X Mn \ ^- . (F4) 

The previous is a typical value, the phase being fixed by 



the microscopic state of the Mn ensemble. Evaluating 
for our parameters we get Xmh ~ 7 x f0~ 3 , typically an 
order of magnitude smaller admixture as that consider 
in the main text, so we can again neglect this admixture 
source. 
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